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Abstract 

Universal Lax pairs of the root type with spectral parameter and independent cou- 
pling constants for twisted non-simply laced Calogero-Moser models are constructed. 
Together with the Lax pairs for the simply laced models and untwisted non-simply 
laced models presented in two previous papers, this completes the derivation of univer- 
sal Lax pairs for all of the Calogero-Moser models based on root systems. As for the 
twisted models based on B n , C n and BC n root systems, a new type of potential term 
with independent coupling constants can be added without destroying integrability. 
They are called extended twisted models. All of the Lax pairs for the twisted models 
presented here are new, except for the one for the F& model based on the short roots. 
The Lax pairs for the twisted G2 model have some novel features. Derivation of vari- 
ous functions, twisted and untwisted, appearing in the Lax pairs for elliptic potentials 
with the spectral parameter is provided. The origin of the spectral parameter is also 
naturally explained. The Lax pairs with spectral parameter, twisted and untwisted, 
for the hyperbolic, the trigonometric and the rational potential models are obtained as 
degenerate limits of those for the elliptic potential models. 



1 Introduction 



This is the third paper in a series devoted to the construction of universal Lax pairs for 
Calogero-Moser models based on root systems. In the first paper paper [Q (this paper will 
be referred to as I), a new and universal formulation of Lax pairs of Calogero-Moser models 
based on simply laced root systems was presented. In the second paper || (this paper will 
be referred to as II), the universal Lax pairs for untwisted non-simply laced models with 
independent coupling constants for the long and short roots were constructed. The twisted 
non-simply laced models were derived from the simply laced models by folding with respect 
to a discrete symmetry of the models. This paper has two objectives: firstly to complete the 
construction of universal Lax pairs of all Calogero-Moser models based on root systems by 
constructing those for the twisted non-simply laced models with independent coupling con- 
stants and secondly to elucidate various functions appearing in the Lax pairs with spectral 
parameter Q, for the elliptic potential. The Calogero-Moser models are a collection 

of completely integrable one-dimensional dynamical systems characterised by root systems 
and a choice of four long-range interaction potentials: (i) 1/L 2 , (ii) 1/ sin 2 L, (ii) 1/ sinh 2 L 
and (iv) p{L), (Weierstrass function) in which L is the inter-particle "distance". The hyper- 
bolic, the trigonometric and the rational potential models are obtained as degenerate limits 
of the elliptic potential model. 

For the general background and the motivations of this series of papers and the physi- 
cal applications of the Calogero-Moser models with various potentials to lower dimensional 
physics, ranging from solid state to particle physics and the Seiberg-Witten theory [||, we 
refer to our previous papers I and II and references therein. 

We address two problems in this paper: the construction of universal Lax pairs with 
independent coupling constants for twisted non-simply laced models and a detailed exposition 
of the functions appearing in the Lax pairs for the twisted and untwisted models. Throughout 
this paper we focus on the elliptic potential case with a spectral parameter. The Weierstrass 
function p(u) is characterised by a pair of primitive periods, which is denoted by {2uji, 2uj 3 } 
(following the convention of 0). As shown in paper II, the twisted models have potentials 
p(u) with the above periods for the long roots and potentials 

p(V»)( w ) {n = 3 for the Gi 

model and n = 2 for the other non-simply laced models) for the short roots, which have 
1/n-th period in one direction {2u>i/n, 2u> 3 }. 

Our main result for the universal Lax pairs for twisted non-simply laced models is very 
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simple. The Lax pairs for the twisted models can be obtained from those for the untwisted 
models (given in paper II) by replacing the untwisted functions {x, y, z, Xd, yd, Zd, x t , y t , z t } 
by the twisted functions {x^ 1//n \ y^/ n \ z^^ n \ x^ n \y^ n \z^ n \ x^ n \y^ n \z^ n ^} whose 
explicit forms will be given in section five. All of the Lax pairs for the twisted models 
presented here are new, except for the one for the F 4 model based on the short roots. The 
verification that the Lax equation is equivalent to the canonical equation of motion and the 
consistency of the Lax pairs is essentially the same as for the untwisted models (paper II) 
and the simply laced models (paper I). As for the twisted models based on B n , C n and 
BC n root systems, a new type of potential terms with independent coupling constants can 
be added without destroying integrability. They are called extended twisted models and 
they have three, three and five independent coupling constants for B n , C n and BC n models, 
respectively. Some new features appear in the Lax pairs for the twisted G2 model. 

Our second objective is to give a detailed exposition of the functions, twisted and un- 
twisted, appearing in the Lax pairs. The untwisted functions were used in previous papers 

1 and II without derivation. Here we derive the explicit forms of these functions starting 
from various sum rules which are necessary and sufficient for consistency of the Lax pairs 
for the elliptic potential. The origin of the spectral parameter is also explained. The Lax 
pairs with spectral parameter, twisted and untwisted, for the hyperbolic, the trigonometric 
and the rational potential models are obtained as degenerate limits of those for the elliptic 
potential models. 

This paper is organised as follows. In section two we recapitulate the essential ingredients 
of the Calogero-Moser models with the elliptic potential. Section three is for the detailed 
exposition of the untwisted functions. Section four gives the root type Lax pairs for all of 
the twisted models. Detailed accounts of the twisted sum rules and functions are provided 
in section five. The final section is devoted to a summary and comments. 

2 Calogero-Moser models with elliptic potential 

In order to set the stage and introduce notation let us start with the definition of the 
Calogero-Moser models based on a simply laced root system A, which is associated with 
a semi-simple and simply-laced Lie algebra g of rank r. The roots a,/3, 7, . . . are real r- 
dimensional vectors and are normalised, without loss of generality, to 2: 

A = {a,/3,7, ...}, a e R r , a 2 = a ■ a = 2, Va G A. (2.1) 
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The dynamical variables are canonical coordinates {V} and their canonical conjugate 
momenta {pj} with the Poisson brackets: 

q 1 ,---,Q r , px,---,Pr, W,Pk} = $j,k, W,q k } = {pj,Pk} = 0- (2.2) 

In most cases we denote them by r dimensional vectors q and p Q, 

q = ( q \ . . . 7 q r ) G R r , p = ( Pl , . . .,p r ) G R r , 

so that the scalar products of q and p with the roots a ■ q, p ■ (3, etc. can be defined. The 
Hamiltonian with the elliptic potential is given by 

w= ^ 2 4s p(a ' 9) ' (2 - 3) 

agA 

in which g is a (real) coupling constant. The Weierstrass function p is a doubly periodic 
meromorphic function with a pair of primitive periods {2ui, 2co>3}, ^s(lv 3 /uji) > 0. We adopt 
the convention that the dependence on the primitive periods is suppressed for simplicity: 



u z z — ' 

m, n 

in which Vt m ^ n is a period 



1 



(v - ) 2 Q, 2 



(2.4) 



&"in, n — 2mui + 2nu)3 

and Yl ' denotes the summation over all integers, positive, negative and zero, excluding 
m = n = 0. Obviously the choice of the primitive periods is not unique. Any modular 
(SL(2, Z)) transformation of the pair of primitive periods is again a pair of primitive periods: 

^1 A ( a \\ a 12 \ ( \ rv ■ ■ o i 

7 = , <iy 6 Z, z,j = 1,2, a n a 2 2 - a 12 a 21 = 1. 

^3 / V 0,21 a 22 / \ U 3 / 

Another characterisation of Weierstrass function is through the invariants: 

p(u) = p(u\2ujx,2uj 3 ) = p(u\{g 2l g 3 }), 

in which 

S 2 = 60]T'fi m 4 n , ^ 3 = 140^'fi m 6 n . (2.5) 



1 For A r models, it is customary to introduce one more degree of freedom, q r+1 and p r +i and embed all 
of the roots in R r+1 . 
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As is well known, the Weierstrass function p(z) satisfies the differential equation 



A[p{u)f 
4(p(u) - 



- 92p(u) - g 3 

e i)(p(«) - e 2 )(p(u) - e 3 ), 



(2.6) 



in which 



ei = p{u t ), e 2 



pM, 



e 3 = p(u 3 ), u 2 



(^i + uj 3 ). 



(2.7) 



It should be noted that the above Hamiltonian is real (hermitian) for real dynamical 
variables p, q and coupling constant g, provided 



One can always choose a pair of primitive periods in this way when the two invariants g 2 , g 3 
are real and the discriminant 



3 Functions in the Lax pairs for untwisted models 

Let us now discuss the functions appearing in the root type Lax pairs for the simply laced 
models and untwisted non-simply laced models. These are the cases treated in Paper I jl| and 
II ||. The elliptic potential case with spectral parameter will be discussed in greater detail 
and the functions appearing in the other cases, the rational, trigonometric and hyperbolic 
cases will be obtained simply as the degeneration of those in the elliptic case. The origin 
of the spectral parameter for the elliptic case is naturally explained in our approach. It also 
leads, as a direct consequence of the degeneration, to a not so widely recognised fact that the 
spectral parameter can be introduced for the rational, trigonometric and hyperbolic cases, 
as well. 

As shown in papers I and II, the Lax pair for the untwisted Calogero-Moser models 
contains functions {x(u),y(u), z(u)}, {xd(u),yd(u), z d {u)} and {x t (u) , y t (u) , Zt(u)} (only in 
the untwisted G 2 model Lax pair based on long roots). They are required to satisfy the first 
sum rule: 



uj\ : real and oos : pure imaginary. 




Dis = gl - 27 'g\ 



(2.9) 



is positive. 



x'{u)x{v) — x'{v)x{u) = x{u + v)[p(v) — p(u)], u,d eC. 



(3.1) 
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The same equations should hold for x d and x t} too. They all have a simple pole at u = 
with unit residue: 

lim ux(u) = \im uxd(u) = \m\ux t {u) = 1. (3.2) 

u^O u—>0 u^Q 

The functions x and Xd are required to satisfy the second sum rule (see (1.2.29)): 

x(u - v ) [p(v ) - p{u)\ + 2 [x d (u)y(-u- v) - y(u + v) x d (-v)] 

+ x(u + v) yd(-v) - Vd{u) x{-u - v) = 0, (3.3) 

and x, Xd and Xt are required to satisfy the third sum rule (see (II. 4. 72)): 

x{3u — 3v)[p{2u — v ) — — 2v)] — x(3f ) — 2v) + yt(2w — u) x(—3u) 
-2x d {3u) y t {—u -v) + 2y t (u + v) x d {-3v) - 3x t (2u - v) y{-3u) + 3y{3v) x t {u - 2v) 
-3x t {u + v) y d {-3v) + 3y d {3u) x t (-u - v) = 0, (3.4) 

in which 

y(u) =z'(iO, j/d(w) = x' d (u), y t (u) = x' t {u). (3.5) 
The z functions are defined as the products of x functions: 

z(u) =x(u)x(-u), z d {u) = x d {u)x d {-u), z t (u) = x t (u)x t (-u). (3.6) 

It should be remarked that the set of solutions {x(u),Xd(u),Xt(u)} to these sum rules has 
symmetry transformations or a kind of 'gauge freedom'. If {x(u), x d {u), Xt(u)} satisfies the 
first, the second and the third sum rules, then 

{x(u) = x(u)e b \ x d (u) = x d (u)e 2bu , x t (u) = x t (u)e 3bu } (3.7) 

also satisfies the same sum rules. Here b is an arbitrary u- independent constant, which can 
depend on the spectral parameter £. The functions z(u), Zd(u), z t (u) are invariant under 
these transformations. 

It is interesting to note that the first ( p.l|) , second ( [3. 3D and third (|3.4| ) sum rules corre- 
spond to the three different types of rank 2 root systems, A 2 , S 2 (C 2 ) and G 2 , respectively. 
The rank 2 root systems are spanned by the roots (3 — k and k — 7 in which k are non-trivial 
intermediate roots in the off-diagonal parts of the commutator: 

[L, M]p tl = ^ {Lp jR M Kil — Mp tK L Kn ) : P,j,keA. (3.8) 

6 



In fact, these three sum rules are obtained as the necessary and sufficient conditions for the 
consistency of the root type Lax pairs in the two-dimensional subspace spanned by (3 — k 
and k — 7. 

Since the first sum rule (|3.1|) is the most fundamental, let us discuss its general solution. 
By taking the limit v — > —u in (|3.1|) and using the residue condition ( |3.2| ) we obtain 



— {x(u)x(— u)} 

Lb Hi 



dp{u) 
du 



From this we know that 



x(u)x(— u) = —p(u) + const. 
We can always choose this constant to be 



(3.9) 



(3.10) 



(3.11) 



and obtain a factorised form ||: 

x{u)x{-u) = -p(u) + p(f ) 



(3.12) 



a{i)a{u)a{i)a{uY 

Here the Weierstrass sigma function o~(u) is defined from p(u) via the Weierstrass zeta 
function ((u) as: 

p(u) = — C(u), ((u) = d log a(u)/du = a'(u)/a(u), 

C(«) = C(^|2^i,2^ 3 ) = - + E ,r ' ' 



+ 



+ 



<* ^ L rn,n ^ L m,n iL m,n 



a(u) = o-(u|2^i,2^ 3 ) = u JJ' (l - — — ) 

\ » 'm, n J 



exp 



u u 
+ 



Q ' 20? 



(3.13) 



in which Yl ' denotes the product over all integers, positive, negative and zero, excluding 
m = n = 0. This is the simplest explanation of why the spectral parameter £ appears in the 
theory. It is elementary to see that 



x(u) = x (u,£), x o( u >€) 



(3.14) 



a{i)a{u) 

satisfies the first sum rule (|3.1| ) for an arbitrary constant £. By the following monodromy 
property of x (u, £): 



x (u + 2(j j ,£) = e 2r >^xo(u,t), x (u,Z + 2u j ) = e 2r » u x (u,£), Vj = (fa), J = 1,2,3 



we find that 

[x' (u, £)x (y, - x' (v, £)x (u, 0] /x Q (u + v, £) 

is a doubly periodic meromorphic function of u, v and £. Meromorphic functions are com- 
pletely determined by their zeros, poles and their residues. It is trivial to check the zeros, 
poles and residues coincide with those of p(v) — p(u). The other sum rules can be verified 
in a similar way. 

Next we analyse possible poles and zeros of x(u). From ( |3.1| ) it is clear that x{u) cannot 
have a pole at u 7^ (modulo periods), because the r.h.s. has no singularity there. From 
the factorisation ( |3.12| ), it also follows that x(u) cannot have a zero at u 7^ {0,±£}. This 
is because an additional zero in x(u) must be canceled by a pole in x(—u) which is not 
allowed as above. By the same argument, we find that x(u) can have no other isolated 
singularities (essential singularities or branch points) at finite u. Thus we deduce that x{u) 
can be expressed as 

x(u) = x (u)e 9( > u \ x (u) = —7— — — , g(—u) = —g( u ) '■ entire function. (3.15) 

o-(6M«) 

By substituting the above form in Q3.1|), we obtain 

x (u)x (v) (g'(u) -g'iv)) 
= x (u + v) (e 9{u+v) - g{u) ~ 9{v) - 1) - p{u)) . (3.16) 

By comparing both sides at u = 00, we find the asymptotic behaviour can only be matched 
for 

I bu b : const. (3-17) 
It should be mentioned that Xq(u) is of exponential growth, i.e., Xq{u) ~ e cu for some constant 
c. The constant b can depend on £. Thus we arrive at the conclusion that the general solution 
of the first sum rule (|3.1| ) is given by 

x{u) = a f c ~f e 6 ", e, b: const. (3.18) 

A similar but less direct analysis of solutions of a functional equation which is a certain 
generalisation of the first sum rule has been given in [|TUJ. 

Since the functions Xd and x t must be of the above general form ( |3.18| ), it is easy to see 



that the general solutions of the first, second and third sum rules are given by (see (II. 2. 20), 
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(II.2.21), (II.4.74)): 

a{^)a(u) a{2^)a(u) <j{3£,)a(u) 



This solution also gives the same z, z d , z t fl3.6|) up to constants: 

z(u) - z d (u) = p(0 - p(20, z(u) - z t (u) = p(0 - p(3£). (3.20) 
It is easy to derive the following relation 

x(2u)x d (-u) + x(-2u)x d (u) = -p{u) + p{£) (3-21) 



from the second sum rule ( |3.3|) . One only has to follow the same path as the factorisation 
identity obtained from the first sum rule. That is, to take the limit v — > u in (|3.3|) 



(x{2u)x d {-u) + x(-2w)x rf (w)) = -— p(u) 
au du 



and integrate it. The above equation and (|3.21|) are responsible for the "renormalisation" of 



the short root coupling constant in the untwisted BC n model (II. 4. 85). In a similar way one 



obtains the following relation from the third sum rule (|3.4|) : 



x(3u)x t (-u) + x(—3u)x t (u) + x d {3u)x t (-2u) + x d {— 3u)x t (2u) = —p(u) + p(£). (3.22) 

At the end of this subsection, let us remark that one can rewrite the first sum rule in an 
equivalent form 

x'{u)x{v) — x'(v)x(u) = x(u + v)[x(u)x(—u) — x(v)x(—v)], (3.23) 
which can be read as an addition theorem: 

x{u + v) = (x'(u)x(v) — x'(v)x(u))/[x(u)x(— u) — x(v)x(—v)]. (3.24) 

3.1 Degenerate cases 

It is well known that the degenerates cases of Weierstrass functions occur when one or both 
of the periods become infinite. The degenerate cases occur if two or all three of e\,e2,e3 
(277) coincide ||. There are three cases, corresponding to the hyperbolic, trigonometric 



and rational potential case, respectively. The degenerate limit of the functions {x, x d , Xt} in 
( 3.19|) gives the corresponding Lax pair with spectral parameter, which is a new result. 
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3.1.1 Real period infinite: Hyperbolic potential 

In this case || 

e\ = e 2 = a e R + , e 3 = —2a, g 2 = 12a 2 , g 3 = —8a 3 , u)\ = oo, cj 3 = —= (3.25) 

yl2a 

and 

p(u\{12a 2 , — 8a 3 }) = a + k 2 s > k = \/3a, 

sinh ku 

C(u|{12a 2 , — 8a 3 }) = — au + /tcoth ku, 

a(u|{12a 2 ,-8a 3 }) = sinh[ ku] exp[— au /2]/k. (3.26) 

The Calogero-Moser model has a hyperbolic potential. The functions appearing in the Lax 
pair for the elliptic potential with spectral parameter £ 

xt(u,0 = =x(u,30, (3.27) 

reduce to 

x(u,£) = k (coth ku — coth k£) e a ^ u , x d (u,£) = k (coth/tu - coth2/«£) e 2a?u , 
x t (n,0 = k (coth ku- coth 3k£) e 3<n . (3.28) 



By utilising the "gauge freedom" ( |3.7|) , the exponential part of the above functions can be 
removed to obtain simple spectral parameter dependent functions: 

x{u, £) = k (coth ku — coth £) = ft(coth km — coth 2k£), 

a? t (w, £) = ^(cothKM — coth3«;£). (3.29) 

If we take the limit — > ±oo, we obtain the spectral parameter independent functions 

x(u) = Xd(u) = x t (u) = a (coth aw =F 1), a : const, (3.30) 

which are essentially the same as the spectral parameter independent functions given in 
previous papers (1.2.7), (11.2.13), (II.4.73). 
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3.1.2 Imaginary period infinite: Trigonometric potential 

In this case 

ei = 2a E R+, e 2 = e 3 = —a, g 2 = 12a 2 , g% = 8a 3 , uj\ = , co> 3 = zoo (3.31) 

Vl2a 

and 



p(«|{12a 2 , 8a 3 }) = — a + k 2 ^-^ — , k — v3a, 



1 

sin 2 km 



C(w|{12a 2 , 8a 3 }) = au + k cot ku, 

a{u\{12a 2 ,8a 3 }) = sin[«w] exp[aM /2]/k. (3.32) 

The corresponding Calogero-Moser model has a trigonometric potential. The functions ap- 
pearing in the Lax pair for the elliptic potential with spectral parameter £ (|3.27|) reduce 
to 

x(u, £) = k (cot km — cot k£) e _ct ^ n , ^(u, £) = « (cot km — cot 2k£) e _2a ^", 
ar t (u,f) = k (cot km - cot 3k£) e" 3<n . (3.33) 

As above the the exponential part of the above functions can be removed by using the "gauge 
freedom" of 



x(u,£) = k(coI km — cot £d(M, £) = K(cot km — cot 2k£), 

Xt(w, £) = /t(cot km — cot 3k£). (3.34) 

If we take the limit k£ — > ±200, we obtain the spectral parameter independent functions 

a; (it) = = ^(m) = a (cot om ± i), a: const, (3.35) 

which are essentially the same as the spectral parameter independent functions given in 
previous papers (1.2.6), (II.2.12), (II.4.73). 

3.1.3 Both periods infinite: Rational potential 

This case is very simple: 

ei = e 2 = e 3 = 0, g 2 = g 3 = 0, u x = -iuj 3 = 00 (3.36) 

and 

p(u\{0, 0}) = m- 2 , C(«|{0, 0}) = u-\ a(u\{0, 0}) = u. (3.37) 
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In this case the Calogero-Moser model has the rational potential and the functions appearing 



in the Lax pair for the elliptic potential with spectral parameter £ (|3.27|) reduce to 



x{u, = - - ~, x d (u, £) = !-—, x t (u, £) = - - J- (3.38) 
M 4 tt 2£ it 3£ 

Thus the Lax pair for the rational potential can also have a spectral parameter, although 
the dependence on it is very simple. By taking the limit as the spectral parameter goes to 
infinity, we recover the well known spectral parameter independent functions 

x(u) = x d (u) = x t (u) = -, (3.39) 

u 

which have been used in our previous papers (1.2.5), (II. 2. 11), (II. 4. 73). 
In the above three degenerate cases, the discriminant ( |2.9| ) is zero. 



4 Universal Lax pairs for Calogero-Moser models based 
on twisted non-simply laced algebras 

In this section we present universal Lax pairs for twisted Calogero-Moser models with in- 
dependent coupling constants. In section 3 of paper II, the twisted Calogero-Moser models 
are derived from the simply laced models with the elliptic potential by folding with respect 
to a discrete symmetry originating from the automorphism of the extended Dynkin dia- 
grams combined with the periodicity of the potential. The integrability of these reduced 
models is inherited from the original simply laced models. In these reduced models, how- 
ever, the coupling constants for the long and short root potentials have a fixed ratio, since 
the simply laced models have only one coupling. In order to exhibit the full symmetry 
of the twisted non-simply laced models, we construct root type Lax pairs with indepen- 
dent coupling constants. As in the untwisted non-simply laced models, there are two kinds 
of root type Lax pairs for the twisted non-simply laced models, the one based on long 
roots and the other short roots. Both are relatively straightforward generalisations of the 
root type Lax pairs for untwisted non-simply laced systems. Our main result for the uni- 
versal Lax pairs for twisted non-simply laced models is very simple. The Lax pairs for 
the twisted models can be obtained from those for the untwisted models (given in paper 
II) by replacing the untwisted functions {x, y,z, x d ,yd,Zd, x t ,y t ,z t } by the twisted func- 
tions {x^^ n \ y^^ n \ z^ n \ x d 1 , y d 1 , z d l ' n , x^ n \y^ n \z^ n ^} whose explicit forms will 
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be given in section five. Among the twisted models, C n , B n and BC n models are differ- 
ent from the rest, the _F 4 and G2 models. In the former models, the long roots of C n 
(±2ej,j = 1, . . . , n), the short roots of B n (±ej,j = 1, . . . , n) and the long and short roots 
of BC n (±2ej, ±ej, j = 1, . . . , n) share one special property, that is, they are orthogonal to 
each other: 

a -13 = 0, a ^{±(3, ±^{3, ±2(3}. (4.1) 

This fact enables us to introduce another independent coupling constant in twisted C n and 
B n models and two more independent couplings in the BC n model. Let us call these models 
with extra independent coupling constant (s) the extended twisted models. 



Except for the Lax pair for the twisted F4 model based on short roots || treated in |4.3.1 
all the Lax pairs reported in this section are new. The verification that the Lax equation 
is equivalent to the canonical equation of motion and the consistency of the Lax pairs is 
essentially the same as in the untwisted models (paper II) and the simply laced models 
(paper I). 

4.1 Extended twisted C n model 

The set of C n roots consists of two parts, long roots and short roots: 

Ac„ = A,UA s , (4.2) 

in which the roots are conveniently expressed in terms of an orthonormal basis of R n : 

A, = {H,T,0,...,} = {±2 ej : j = l,...,n}, 2n roots, (4.3) 

A s = {a,/3,j,...,} = {±ej±e k : j ^ k = 1, . . . , n}, 2n(n - 1) roots. (4.4) 



The Hamiltonian of untwisted C n model is given by 

HZ wisted = \v 2 + I £ P(a ■ q) + f ^ S • 9), (4-5) 



2 1 2 ^ y 11 4 ^ v 

in which gi and g s are the two coupling constants for the long roots and short roots. The 
twisted model is obtained by multiplying one of the primitive periods of the potential for 
the short roots by a factor of one half (II. 3. 21): 
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In this paper we choose {2ui, 2u>s} — > {u>\, 2u>s}, that is the potential for the short roots is 

p {1/2) (u) = p(u\uj u 2uj 3 ) = p(u) + p(u + uji) - p(u x ). (4.7) 

The second equality is known as Landen's transformation. The other case {2uji,2uj 3 } — > 
{2u>i, u 3 } can be treated in a similar way. The extended twisted C n model is defined, roughly 
speaking, by a sum of the untwisted Hamiltonian at a half-period and the twisted Hamilto- 
nian: 

n_i ex— twisted njuntwisted (1/2) , njtwisted 

rt r = rtr> + rln 

= V + f E ^ (1/2) ( a • + 1 E • ?) + &p {m & ■ i)] - ( 4 - 8 ) 

in which the "renormalised" long root coupling constant is defined by 

~9l=9h{9h + ^9h) (4-9) 



from the parameters in the Lax pairs Q4.14 ) and ( |4 .19 ). The potential for the long root has 



two terms, the Weierstrass function with periods {2a>i, 2uj^} and the one with the half-period 
{o;i,2o;3}. Obviously the above Hamiltonian (|4.8| ) is integrable for gf = or gf = 0. The 
orthogonality of the long roots fl4.1|) allows these two potentials to co-exist without breaking 
the integrability. 

The Lax pairs for the twisted model are simply obtained from those of the untwisted 
model by changing the functions x, Xd and y, yd corresponding to short roots to new ones 
having a half-period. 

4.1.1 Root type Lax pair for extended twisted C n model based on short roots 

The Lax pair is given in terms of short roots. It is determined by the root difference pattern 
of the short roots: 

short root 

,, . . i 2 x short root , , 1 _s 

C n : short root — short root = < , (4.1U) 

' long root 

non-root 

The matrix elements of L s and M s are labeled by indices 7 etc.: 

L.{q,p,£) = p-H + X + X d + X L , 
M s (q,0 = D + D L + Y + Y d + X L , (4.11) 
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Here X and Y correspond to the part of "short root — short root = short root" of ( |4.10|) : 



X = ig s ^ (V2) (« • q,OE(a), Y = tg s J2 2/ (1/2) (« • 9, E(a) M = V 7 ,a, 

(4.12) 

and Xd and Yd correspond to "short root — short root = 2x short root" of (|4.10|) : 

X d = 2ig s ^ x d ( a ■ qiO E d(a), Y d = ig s ^ Vd^i® ' q^) E d (a), E d {a) Pl = 8^ M - 

(4.13) 

The terms in L s (M s ), X^ (Yl) corresponding to "short root — short root = long root" of 
( |4 .17 ) have two terms corresponding to two independent coupling constants for the long root 



potentials: 

X l = ^MS-g,0+5i/ /2) (2'g,0]^), 
HeA; 

Yl = i y £[9hy(Z-q,0+9iJ l/2) (Z-q,(i)]E(E), 
hga ; 

E(E) Pl = V 7)3 . (4-14) 
The diagonal parts of L s and M s are given by 

H fh = /36 fi „, D fh = 6 fi „D fl , Dp = -zg s (z^ 2 \f3-q)+ £ ^ 1/2) ^-q) 

\ KeA s , k-/3=1 

(4.15) 

and 

(D L ) M = 5p„(D L ) p , {D L ) p = -i Yl ■ Q) + 9i 2 z {l/2) (T ■ q)] . (4.16) 

TeA;, /3-T=2 

The new functions x^ 1//2 - ) , f/ 1 / 2 ) = a^ 1 / 2 )', z^l 2 ^ and x^ 2 \y^ 2 ^ = x^ 2 ^ must satisfy certain 
sum rules which are generalisations of the second sum rule. These sum rules are listed in 
section |5] together with the explicit forms of the new functions. It is easy to verify 

Tr(L 2 ) = 8(n - i)n e c x ~ twisted , 



in which a factorisation identity ( 5.23Q is responsible for the renormalisation of the long root 
coupling. 
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4.1.2 Root type Lax pair for extended twisted C n model based on long roots A/ 



The Lax pair is given in terms of long roots. It is determined by the root difference pattern 
of the long roots: 



2 x long root 

C n : long root — long root = { 2 x short root 

non-root 



(4.17) 



The matrix elements of Ll and Ml are labeled by indices T, Q etc.: 

L L (q,P,0 = p-H + X d + X s , 
M L (q,0 = D + Y d + Y s . (4.18) 

Here X d and Y d correspond to the part of "long root — long root = 2x long root" of (|4.17|) : 

E d (E), 



EJE) 



seA ; 

Y d = i ^ [9hVd{E ■ q, + 9i 2 y d /2) (2 • q, £) 
seA, 

E d (^)m = ^T-f2,2H, 
and X s and Y s correspond to "long root — long root = 2x short root" of ( f4.17| ): 



(4.19) 



X s = 2ig s ^ x [ 1 ( a -q,£)E d (a), Y s = ig s ^ Vd ( a 'Q>€) E d(oi), E d (a) rn = £ T -<>.2„- 



aeA s 



(4.20) 



The diagonal parts of Ll and Ml are given by 
D r = 



-i I g h z d (T ■ q,£)+ g h z d 1/2) (T -q,£)+g s 

V reeA a , k-T=2 



.a/2), 



«■?) ■ (4.21) 



The functions x^^ 2 \ t/ 1 / 2 ) = a^ 1 / 2 )', z 1 ^ 1 / 2 ^ and x d l ^ 2 \y d 1 ^ 2 ^ = x^ 2 ^ are the same as those 
appearing in the Lax pair based on short roots. Needless to say, the consistency of the root 
type Lax pairs (|4.11|) and ( |4.18| ) does not depend on the explicit representation of the roots 
in terms of the orthonormal basis ( |4.3| , |4.4| ). This remark applies to the other models as well. 



4.2 Extended twisted B n model 

The set of B n roots consists of two parts, long roots and short roots: 

A Bn = A t U A a , 



(4.22) 
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in which the roots are conveniently expressed in terms of an orthonormal basis of R n : 

A, = {a,P,j, ...,} = {±ej ± e k : j ^ k = 1, . . . , n}, 2n(n - 1) roots, 

A s = {A, [A, u, ...,} = {±ej : j — l,...,n}, 2n roots. (4.23) 

The Hamiltonian of the untwisted B n model is given by 

HZ tmsted = \p 2 + | E P(a ■ q) + g] £ p{\ ■ q). (4.24) 

As in the C n case, the twisted model is obtained by multiplying one of the primitive periods 
of the potential for the short roots by a factor of one half (II. 3. 33): 

= \v 2 + f E P(a ■ q) + gl Y, ^ (V2) (A • q). (4.25) 

agA( AgA s 

The extended twisted B n model is defined, roughly speaking, as a sum of the untwisted 
Hamiltonian and the twisted Hamiltonian: 

ni ex— twisted njuntwisted i n/twisted 

tip, = n R + tin 

= \v 2 + f E P(« ■ ?) + E ^^( A • ?) + 9 2 s J 1/2 \\ ■ q)} , (4.26) 

oeA; AeA s 

in which the renormalised short root coupling constant is defined by 

~g 2 1 =g sl (9s 1 +2g S2 ) (4.27) 
from the parameters in the Lax pairs ( 4.31]) and (|4.38 ). 



4.2.1 Root type Lax pair for extended twisted B n model based on short roots 

The Lax pair is given in terms of the short roots. It is determined by the root difference 
pattern of the short roots: 

{long root 
2 x short root (4.28) 
non-root 

The matrix elements of L s and M s are labeled by indices /i, v etc.: 

L s (q,P,0 = p-H + X + X d , 
M s (q,0 = D + Y + Y d . (4.29) 
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Here X and Y correspond to the part of "short root — short root = long root" of (|4.28|) : 

X = ig t ^2 x ( a ' E ( a )> Y = i 9i^2 y ( a ' g ' €) E (. a )i E(a)n V = <V-„ jQi . (4.30) 

Thus they have the same form as in the untwisted model. The difference occurs in X d and 
Y d , which correspond to "short root — short root = 2x short root" of ( f4.28|) : 

X d =2iJ2 [9.iXd(X-q,$+gs a x$ /2) (X'q,£)\E d (\), 
AeA s 



AeA s 



E d (X) 



5. 



(4.31) 



Both full and half-period functions appear in these matrix elements. The diagonal parts of 
L s and M s are given by 



H, 



D, 



\ 7-U=l / 



The functions z^ 1 ! 2 ^ and 



(1/2) (1/2) _ (1/2)' 



are the same as those that appear in the Lax 



pairs for the twisted C n models. It is easy to verify 

Tr(L 2 s ) = 4<H%- twisted , 



(4.33) 



in which the extended twisted B n Hamiltonian is given above ( 4.26 ) 



4.2.2 Root type Lax pair for extended twisted B n model based on long roots A; 

The Lax pair is given in terms of the long roots. The matrix elements of L\ and M\ are 
labeled by indices a, (3 etc. From the root difference pattern of the long roots 



B r . 



we obtain the Lax pair: 



long root — long root 



long root 
2 x long root 
2 x short root 
non-root 



(4.34) 



Li(q,P,0 = P-H + X + X d + X s , 
M,(g,0 = D + Ds + Y + Y d + Y s . 



(4.35) 
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Here X and Y correspond to the part of "long root — long root = long root" of (|4.34 ): 

X = igi^2 x ^ a ' €) E ( a )> Y = i 9i^2 y ( a ' g ' E ( a )Pi = S P-y,a, (4-36) 

and X d and Y d correspond to "long root — long root = 2x long root" of (|4.34j) : 



X d = 2ig t ^ x d (a-q,£)E d (a), Y d = ig x ^ Vd{oc-q^)E d (a), E d (a) Pl = 5/3^ 7)2a . (4.37) 

An additional term in Li (Mi), X s (Y s ) corresponds to "long root — long root = 2x short 
root" of ( fOg) : 



x s = 2tJ2 [g sl Xd^-q,0 + 9s 2 x d 1/2) ^-Q^)}E d (X), 

AeA s 
aga s 

E d (X)i3j = <5/3- 7 ,2A- 

The diagonal parts of L/ and M/ are given by 

#s 7 = /% 7 , D/3 7 = 5/3,7^/3, £>/3 = -igi lz(p-q)+ Z ( K ' q "> ) 

V kGA;, k-/3=1 / 



(4.38) 



(4.39) 



and 



(Ds) fh = 6 fltl (D8)f i , (Ds)p = -i J2 [9sM^-Q,0+9sA 1/2) ^-^0\ • (4-40) 

AeA s , /3-A=l 

The functions z^ 1 / 2 ) and x^ 2 ^ , y^ 2 "* = x^ 2 ^ are the same as those that appear in the Lax 
pairs for the twisted C n models. It is easy to verify that 



Tr(Ll) = 8(n - 



ex— twisted 



(4.41) 



in which the extended twisted B n Hamiltonian is given above (|4.26| ) 



4.3 Twisted F 4 model 

The set of _F 4 roots consists of two parts, long and short roots: 



A F4 = A t U A s 



(4.42) 
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in which the roots are conveniently expressed in terms of an orthonormal basis of R : 



A s = = {±e j ,-(±e 1 ± e 2 ± e 3 ± e 4 ) : j = 1, ...,4}, 24 roots.(4.43) 



A/ = {a, (3, 7, . . . , } = {±ej ± e k : j ^ A; = 1, . . . , 4}, 24 roots, 

1 

2 l 

The set of long roots has the same structure as the D4 roots and the set of short roots has 
the same structure as the union of -D4 vector, spinor and anti-spinor weights. The untwisted 
F4 Hamiltonian is given by 

Kg**** = l p 2 + I p(a ■ q) + g 2 s £ p(X ■ q). (AAA) 

aeAt XeA s 

It should be noted that this has the same general structure as the Hamiltonian of the 
untwisted B n theory Q4.24 ). As in the C n and B n cases, the twisted model is obtained 



by multiplying one of the primitive periods of the potential for the short roots by a factor 
of one half (II.3.41): 

KT ted = \ V 2 + f £ p{a ■ q) + gl P (V2) (A • g). (4.45) 
qgA; AeA s 

4.3.1 Root type Lax pair for twisted F 4 model based on short roots A s 

The Lax pair is given in terms of the short roots. It is determined by the root difference 
pattern of the short roots: 

{long root 
s ort root (4.46) 
2 x short root 
non-root 

The matrix elements of L s and M s are labeled by indices /x, v etc.: 

L a {q,P,0 = vH + X + X d + X h 
M s (q,0 = D + D l + Y + Y d + Y l . (4.47) 

Here X and Y correspond to the part of "short root — short root = short root" of ( |4.46|) : 



X = tg s J2x il/2 \\-q,OE(\), Y = ig s ^ y {1/2) (X-q, £)E(\), E{\), v = 5,. v ^ (4.48) 
AeA s AeA 3 

and Xd and Yd correspond to "short root — short root = 2x short root" of ( |4.46| ): 

X d = 2ig s 4 V2) (A ' q,C)E d (X), Y d = ig s ^ y d 1/2 \\ ■ q,£)E d (\), E d {\), v = V„ |2A , 
AeAs AeAs 

(4.49) 
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and these functions have half-period. The additional terms X\ and Y\ correspond to "short 
root — short root = long root" of ( |4.46 ) and they have the same form as in the untwisted 



model: 

Xi = igi ^2 x ( a ' 9> €) E ( a )> Y i = ^ y ( a ' 9 > €) E ( a )> E ( a )^ = (4.50) 
The diagonal parts of L s and M s are given by 

H^u = fJiS^, = S^Dp, Dp = -ig s I z^ (ji-q)+ ^ ^ (1/2) (A • q) 

\ XeA s , X-n=l/2 

(4.51) 

and 

(% = ^(% {Di)„ = -igi Yl ( 4 - 52 ) 

CtGA;, a-fi=l 

The functions x^ l ^ 2 \ y^^ 2 \ z^ 1 / 2 ) and x} \y£ are the same as those appear in the Lax 
pairs for the twisted C n models. It is easy to verify that 

Tr(L 2 ) = 12H% isted , (4.53) 

in which the twisted Hamiltonian is given by ( |4.45| ). In different notation the above Lax pair 



( |4.47| )- (f4.52|) was reported in 



4.3.2 Root type Lax pair for untwisted F 4 model based on long roots A; 

The Lax pair is given in terms of long roots. The general structure of this Lax pair is 
essentially the same as that of the B n theory, since the pattern of the long root- long root 

long root 

i i j 2 x long root , , . . 

/ , : long root - long root = { g x ^ ^ (4.54) 

non-root 



is the same as that of B n ( |4.34j ). So we list the general form only without further explanation. 



They are matrices with indices /3, 7 etc.: 

Li(q,P,0 = V H + X + X d + X s , 
Af,(g,0 = D + Ds + Y + Y d + Y s . (4.55) 

X = igi^2 x(a ■ q, £)E(a), Y = igi ^ y(a ■ q, £)E(a), E(a) Pl = fy-^a- (4.56) 
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X d = 2%gi ^ x d (a-q,^)E d (a), Y d = ig x ^ yd(a-q,£)E d (a), E d (a) Pl = ^- 7 , 2a - (4-57) 
X s = 2ig s 4 V2) (A " q,Z)E d (\), Y s = ig s £ y d 1/2 \\ ■ q,Z)E d (\), E d {\) Pl = ^_ 7>2 , 



aga s aga s 



The diagonal parts of L x and Mi are given by 



H/3y = f35p }1 , D Pl = dfrjDp, Dp = -ig t I z(j3 ■ q) + ^ z(k ■ q) J 

\ kGA;, k-/3=1 / 



(4.58) 



(4.59) 



and 

(Da) fh = 5 fin {p8)f ) , {Ds)p = -ig s ]T z^ 2 \\-q). (4.60) 

AGA S , /3-A=l 

The functions z^ 1 / 2 ^ and x^^yU^^ are the same as those that appear in the Lax pairs for 
the twisted C n models. It is easy to verify that 

Tr{Lf) = 2AH% isted , (4.61) 

in which the twisted _F 4 Hamiltonian is given above (|4.45|) . 

4.4 Twisted G 2 model 

A Lax pair for the twisted G 2 model has never been constructed in any form. Here we 
present new Lax pairs based on short and long roots for this model. The twisted G 2 model 
has some features different from those discussed above due to the ratio of (long root) 2 / (short 
root) 2 = 3 instead of 2 for the other non-simply laced root systems. The set of G 2 roots 
consists of two parts, long and short roots: 

A G2 = A, U A s , (4.62) 

in which the roots are conveniently expressed in terms of an orthonormal basis of R 2 : 

A* = {a,/?, 7 ,...,} = {±(-3ei + V3e 2 )/2,±(3 ei + V3e 2 )/2,±V3e 2 }, 6 roots, 
A s = {X,n,u,...,} = {±e 1 ,±(-e 1 + V3e 2 )/2,±{e l + y/3e 2 )/2}, 6 roots. (4.63) 

The sets of long and short roots have the same structure as the A 2 roots, scaled and rotated 
by 7r/6. The untwisted G 2 Hamiltonian is given by 

n^ wi sted = \p 2 + jY1 ■ 9) + 9 2 s E p( a • ( 4 - 64 ) 

oga ( aga s 
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The twisting leads to a short root potential with a third-period function (II. 3. 48): 

H% sted = \v 2 + | £ p(a ■ q) + gl ^ (1/3) (A " ?), (4.65) 

aeA ; AeA s 

in which we choose the convention ( |4.7|) 

p d/3)( u ) = p( M | 2wi /3,2^ 3 ) 

= p(w) + p(u + 2^/3) + + 4^/3) - p(2wi/3) - p(4w x /3). (4.66) 

The second equality is Landen's transformation for a third-period. 

4.4.1 Root type Lax pair for twisted G 2 model based on short roots A s 

The Lax pair is given in terms of short roots. The general structure of this Lax pair is 
essentially the same as that of theory, since the pattern of the short root- short root 

long root 

i i ; short root , , „_ N 

Go '■ snort root — snort root = < _ , 4.67 

4 2 x short root v ' 

non-root 

is the same as that of the F4 model (|4.46|) . The difference is that the functions associated 
with the short root potential have a third-period. The matrix elements of L s and M s are 
labeled by indices /x, v etc.: 

L,{q,p,0 = V H + X + X d + X h 
M s (q,0 = D + D l +Y + Y d + Y l , (4.68) 



X = tg s J2x^ 3 \X-q,OE(X), Y = ig s £ yW 3 \\-q,£)E(\), E(X) flu = 5^ x . (4.69) 

AeA s AeA s 

X d = 2ig s 4 1/3) (A • q,Z)E d (\), Y d = ig s £ ^ 1/3) (A • q,£)E d (\), E d (\), v = 5^ U;2X . 
AeAs AeAs 

(4.70) 

Xi = igi Y x ( a ' ?> €)E(a), Y t = ig l ^ vi a ' ?> E (a), E(a)^ v = 5^ u>a . (4.71) 

aeA; cteA; 

The diagonal parts of L s and M s are given by 



= fi6^, = 6^D„ D, = -ig s ( z^{p, ■ q) + ^ z^ 3 \\ ■ q) 

AeAs, X-fi=l/2 



(4.72) 
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and 

(AU = <W(A) M) (Di)p = -igi < a -^- ( 4 - 73 ) 

aeAi, a-fi=3/2 

The new functions x^^ 3 \ y^^ 3 \ 2^ 1//3 - ) and x d ,y d 3 ^ must satisfy certain sum rules which 
are generalisations of the third sum rule. These sum rules are listed in section |5] together 
with the explicit forms of the new functions. It is easy to verify that 

Tr(L 2 s ) = 6H^ sted , (4.74) 

in which the twisted G2 Hamiltonian is given by (|4.65|) . 

4.4.2 Root type Lax pair for untwisted G* 2 model based on long roots A; 

This Lax pair is different from the others because of the 'triple root' term in the pattern of 
the long root- long root: 

long root 

1 2 x long root ,, __ N 

G 2 : long root - long root = < g ^ ^ ^ (4.75) 

non-root 

The matrix elements of Li and Mi are labeled by indices /3, 7 etc.: 

Li(q,p,0 = p-H + X + X d + X t , 
M,(?,£) = D + Dt + Y + Y d + Y t . (4.76) 

The terms X, Y and X d , Yd are the same as before: 

X = igi^ X ( a ' E ( a )> Y = i 9l^2 ' g ' E ( a )lh = 8p-y,a- (4-77) 

X d = 2igi x (i (a-g,0-Erf(a), *d = *0z ^ yd(a-q,^)E d (a), E d (a) M = ^- 7)2a - (4.78) 

The 'triple root' terms X 4 and Yt are associated with the new functions x^ 3 ^ and y^ 1 : 

X t = 3ig s ^! 1/3) ( A • Q,OEt(X), Y t = ig s ^ y J 1/3) (A • q,0E t (\), ^(A)^ = 6 P ^ M . 
AeA s AeA s 

(4.79) 

The diagonal parts of L\ and M; are given by 

Epn = D Pl = 8p n Dp, Dp = -ig x I z(/3 ■ q) + z(k ■ q) , (4.80) 

y «eAj, k-/3=3/2 / 
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and 

{Dt)f h = 8 fl „(Dt)p, {Dt)p = -ig s Yl ^ (1/3) (A-g). (4.81) 

XeAs, /3-A=3/2 

The new functions x^ 1 and must satisfy certain sum rules which are generalisations 
of the third sum rule. These sum rules are listed in section |5| together with the explicit forms 
of the new functions. 

4.5 Extended Twisted BC n root system Lax pair with five inde- 
pendent couplings 

The BC n root system consists of three parts: long, middle and short roots: 

A BCn = A t U A m U A a , (4.82) 

in which the roots are conveniently expressed in terms of an orthonormal basis of R n : 

A* = {S,T,f2, ...,} = {±2ej : j = l,...,n}, 2n roots, (4.83) 

A m = {a,/3,j, ...,} = {±ej ± e k : j, k = 1, . . . , n}, 2n(n - 1) roots, (4.84) 
A s = {X, ...,} = {±ej : j = l,...,n}, 2n roots. (4.85) 



The Hamiltonian of the untwisted BC n model is the C n Hamiltonian (|4.5|) plus the contri- 
bution from the short root potential: 

H unt : iste d = 1^ + A J- P (a ■ q) + I ^ S • 9) + 9" E ^ A ' «)" ( 4 " 86 ) 

«eA m seA; agA s 

The Hamiltonian of the extended twisted BC n model is, roughly speaking, the sum of the 
extended twisted C n Hamiltonian (4J) and the extended twisted B n Hamiltonian Q4.26 ) at 
half-period: 



aeA m 5eA ; 

+ E Kp (1/2) (A-?)+^ 2 P (1/4) (A-g)], (4.87) 



AeA s 

in which renormalised coupling constants are defined by 



<?h = 9h{9h + 2g. 



9 2 Sl = 9n [gs! + ~{9h +9h) + 2 9s 2 J + 7}9h9s2, 

9 2 s 2 = 9s 2 (9s 2 + 2g h ). (U 
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The twisted BC n Hamiltonian derived by reduction in paper II (II. 3. 61) may be obtained by 

9h = 9m = g si = -g S2 = 9h = °- 

Here we consider the Lax pair based on the middle roots only. The pattern of middle 
root- middle root is 

long root 
middle root 

2 x middle root (4.89) 
2 x short root 
non-root 



BC n : 



middle root — middle root 



We can then construct the root type Lax pair for BC n root system: 

L m {q,P,Z) = p-H + X + X d + X L + X s , 
M m (g,0 = D + D L + Y + Y d + Y L + Ds + Y S . 



(4.90) 



The matrix elements of L m and M m are labeled by indices (3, 7 etc. This Lax pair is, roughly 
speaking, the sum of the extended twisted C n Lax pair based on short roots (with renaming 
of the couplings) and the extended twisted B n Lax pair based on long roots (both periods of 
the potentials are halved). Here p ■ H + X + X d + X L (D + D L + Y + Y d + Y L ) is exactly the 
same as the L s (M s ) matrix of extended C n models with three coupling constants based on 
short roots. An additional term in L m (M m ), X s (Y s ) corresponds to "middle root — middle 



root = 2x short root" of ( fL89]) : 

X s = 2iJ2 

AeA s 



g Sl x d ll2) (A • q, + 9s 2 x d m (A " q, E d (X), 



{Ds) Pl = 8p a (Ds)p, (Ds)p 



Y s = l ^[9sJ d /2 \\-q,0+9s 2 y d m (X-q,0\E d W, 
AeA s 



(4.91) 



(m (\-q^)+9sJ d /4 \X-q,C)], (4.92) 



AeA s , /3-A=l 

which is a half-period version of the extended twisted B n model based on long roots ( |4.38| ). 

The new functions x d X ^\ y d 1 \ and their sum rules are given in section 5. It is easy 
to verify that 

Tr(L 2 m ) = 8(n - l)^j", (4.93) 

for which various factorisation relations ( |3.21| ),( [5.22j )-( |oT2^| ), are responsible for the renor- 
malisation of the couplings ( 4.88|) . The verification that the Lax equation is equivalent to 
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the canonical equation of motion is again essentially the same as in the extended twisted C n 
and B n models. 

At the end of this section let us remark on the relation between our extended twisted 
BC n model and Inozemtsev's work ||. His model has a potential, which in our notation is: 

v Inozemtsev = y E p {m (*-<i) 

a€A m 

+ \ E [siV V2) (A • ?) + S 2 V 1/2) (A -q + f ) + S 3 y i/2) (A • q + u 3 ) 
AeA s 

+ ^p (1/2) (A-g+y + c 3 )], (4.94) 
which is equivalent to our Hamiltonian ( |4.87| ) if the coupling constants are related by 



~2 of 2 2\ 2 o 2 2 

9h = %3 - 9a), 9h = 8 94 , 9m = 9 

f 2 2 2 , 2\ ~2 * 

2\9i ~ 92 ~ 9a + 9a), 9 S2 = 2 



~g 2 si = ^(gl-gi-gi + gl), ~gl 2 = ^{gl- £)• (4.95) 



5 Functions in the Lax pairs for twisted models 

As we have seen in the previous section, many different functions enter in the Lax pairs for 
the twisted Calogero-Moser models. They are untwisted functions 

x, y = x', z, x d , y d = x' d , z d 

associated with the potential for the long roots. Another set of functions with half-period 

-(1/2) ..(1/2) _ (1/2)' (1/2) (1/2) (1/2) _ (1/2)' (1/2) 

x , y — d, , 4, , x d , y d — x d 7 *d 

is introduced for potential for short roots having the ratio (long root) 2 : (short root) 2 = 2:1. 
They appear in the twisted C n , B n , F 4 and BC n models. The third set of functions with 
third-period 

T (l/3) „(l/3) _ (1/3)' (1/3) (1/3) (1/3) _ (1/3)' (1/3) 

(1/3) (1/3) _ (1/3)' (1/3) 

x t 7 yt ~ x t 7 *t 
is used in the twisted G2 model in which the ratio (long root) 2 : (short root) 2 is 3 : 1. Finally 
the functions with quarter period 

(1/4) (1/4) _ (1/4)' (1/4) 

x d 7 yd — x d 7 ^d 

appear in the twisted BC n model. 
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5.1 Twisted sum rules 



In this subsection we give the sum rules or the functional relations to be satisfied by these 
functions. The explicit forms of these functions will be given in the next subsection. First 
of all, they all have to satisfy the first sum rule: 

x^ ln) \u)x {l/n \v)-x^ ln) \v)x {l/n \u) = x {l/n \u + v)[p^ /n \v)-p {l/n \u)}, 71 = 2,3, (5.1) 

4 d 1/n) \u)x d 1/n \v)-x d l/n) \v)x d 1/n \u) = x d 1/n \u + v)[p^ n \v) - p {1/n) (u)}, n = 2,3,4, 

(5.2) 

xl 1/3)/ (w)xl 1/3) (v) - x^ ] \v)xf /z) (u) = xJ 1/3) (u + v) [p^ (v) - p (1/3) (u)] , (5.3) 
so that they give the potential with the 1/n-th period by the factorisation formula ( [3.10| ): 





[u) 


= x^ 


(u 


)x^ 


-u) 


= -p (1/n \u) + const, 


n = 2,3, 




» 


= 4 ln) 


(t* 


>4 1/n) 


-u) 


= -p {1/n \u) + const, 


n = 2,3,4, 


4 1/3) 


[u) 


= 4 1/3) 


u) 


4 1/3 \ 


-u) 


= -p {1/3) (u) + const. 





There are three types of twisted second sum rules: 

x(u-v) [p {1/2) (v) - p (1/2) (u)] + 2 \x ( d /2) (u)y(-u-v)-y(u + v) x ( d ,2 \-v) 

+ x(u + v) y^ 2 \—v) — y^ 1 (u) x(—u — v) = 0, 



(5.4) 



x {1/2 \u-v)[p{2v) - p{2u)} + x{2u)y {1/2) {-u-v) - y {l/2 \u + v) x(-2v) 

+ x (1/2) (u + v) y(-2v) - y(2u) s (1/2) (-u - v) = 0, (5.5) 

x J /2) (u - v) \p{2v) - p{2u)} + x d {2u) y ( d /2 \-u - v) - y { d m {u + v) x d {-2v) 

+ x d 1/2 \u + v) y d (-2v) - y d (2u) x [ d 1/2 \-u - v) = 0, (5.6) 



in which (O) is obtained from the untwisted second sum rule ( |3.3j ) by changing the double 



root functions x d ,y d ,z d into x d l ^ 2 \ y d , ■ This is related to the two dimensional root 



system B n and a part of its root difference patterns ( |4.28| ), Q4.34| ), i.e. 



(long or short) root — (long or short) root 



long root 

2 x short root (x 



[x) 

(1/2) > 



(5.7) 



This sum rule is responsible for the consistency of the twisted Lax pairs of the models 
containing the above pattern, see subsections |4.1.1| , |4.2.2| , 4.2.1| , [4.3. 1| , |4.3.2j . The next 
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twisted second sum rule (|5.5|) is related to the two dimensional root system C n and to a part 
of its root difference patterns ( |4.f 0| ), i.e. 

short root — short root = < l° n g ro °t (5.8) 

[ short root (x ( ' >) 

This sum rule is responsible for the consistency of the twisted Lax pairs of the models con- 
taining the above pattern, that is the C n and models based on short roots, see subsections 
[4.1.1| , [4.3. 1| . The last twisted second sum rule ( |5.6| ) is related to the two dimensional root 
system B n {C n ) and to a part of its root difference patterns ( [4.17| ), ( |4.34| ), i.e. 



f 2 x long root (xd) 
long root - long root = j ^ ^ ^ ^ ( ^ 1/2)) (5.9) 

Obviously ( |5.6| ) can be obtained from ( |5.5| ) by changing x, x^ 1//2 - ) (and y,y^^ 2 ^ ) to Xd,x^ 2 ^ 
(and y^y^ 2 ^)- This sum rule is responsible for the consistency of the twisted Lax pairs 
of the models containing the above pattern, that is the C n , B n and _F 4 models based on 
long roots, see subsections [4.1.2| , |4.2.1| , |4 . 3 . 2| . Finally, the twisted functions also satisfy the 
ordinary second sum rules: 

x (1/n) (u - v) [p {1/n) (v) - p (1/n) (u)] + 2 \x d 1/n) (u) y {l/n \-u - v) - y {1/n) (u + v) x d 1/n) (-v) 
+ x {1/n) (u + v)y d 1/n \-v)-y d 1/n) (u)x {1/n) (-u-v) = 0, n = 2,3,4. (5.10) 

This is responsible for the consistency of the twisted Lax pairs of the models containing the 
pattern 

i / .ill \ i /•-,!, \ f short (middle) root (a;* 1 /™)) 

short (middle) root — short (middle) root = ^ ^ t ( _ . Jjn N , (i/n)\ (5.11) 

d 



2 x short (middle) root (x 



This occurs in F 4 (n = 2) and G*2 (n = 3) Lax pair based on short roots and BC n (n = 2) 
Lax pair based on middle roots. Some of the above twisted second sum rules, (|5.4j) , (|5.5| ) 
and (|5.10|) for n = 2 were reported in ||. 

There are two types of twisted third sum rules: 

x(3u - 3v)[p ( - 1/3) (2u -v)- p (1/3) (u - 2v)} - x(3v) y t (1/3) (u - 2v) 
+yl im (2u - v) x(-3u) - 2x d (3u) yl 1/3) (-u - v) + 2yf l *\u + v) x d (-3v) 
-3xJ 1/3) (2m - v) y{-3u) + 3y (3v) a^ 1/3) (u - 2u) - 3^ 1/3) (u + v) jfc(-3r;) 

+ 3y d (3w) xJ 1/3) (-u - u) = 0, (5.12) 
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x (1/3) (w - v)[p{2u -v)- p{u - 2v)\ - x {1/3 \v) y{u - 2v) 

+y{2u - v) x (1/3) {-u) - 2xi 1/3) (u) y{-u - v) + 2y(u + v) x { j /3) (-v) 

-x(2u - v) y (1/3 \-u) + y (1/3 \v) x(u - 2v) - x(u + v) y£ /3 \-v) 

+ y% /3) (u)x(-u-v) = 0, (5.13) 



in which ( 5.12|) is obtained from the untwisted third sum rule (|3.4j) by changing the triple root 
functions x t ,y t ,z t into , yjr , -4 ■ Obviously this is related to the two dimensional 



root system G2 and its long root -long root pattern ( |4.75| ). This is responsible for the 



consistency of twisted G2 Lax pair based on long roots discussed in |4.4.2| . The next twisted 
third sum rule ( |5.13|) plays a role in twisted G2 Lax pair based on short roots discussed in 

m 



5.2 Twisted functions 

Here we give a simple set of functions with spectral parameter satisfying all of the sum rules 
given in section 3, (p|), (p|), (Q and section 5, (|5TT0|) , (|5TT2|) and ( gT3|) : 



= -tIw~T' x d(u,€)=z(u,2£), x t (u,£) = x(u,3£), (5.14) 
X(1/2)(M '° = JjW^^) eXP[e ^ U/2] ^ 1/2) ^0 = ^ (1/2) K20, (5.15) 

= ff (l/3)j3) ff (l/3)( M ) eX Pt( 2 / 3 )^ 2 ^/ 3 )^]' 

xi 1/3) ( M ,0 = * (1/3) K2£), s^(Ti,O = z (1/8) (u,30, (5.16) 

* (1/4) («,o = J^^^ ex p [(ei/4 + ^/2)/2) e«] , 

a£ /4) (u,0 = z (1/4) (u,2£), (5.17) 

in which the superscript (1/n) on the Weierstrass sigma functions denotes that their primitive 
periods are {2a>i/n, 2u> 3 } 

a {1/n \u) = a{u\2uJi/n,2LOz) 

and those without them have the usual periods {2uji, 2u 3 }, i.e., have n — 1. In order to 
show that they actually satisfy the sum rules, the monodromy property of these functions 
is important. The minimum intervals and corresponding monodromies in u and £ for these 
functions are given in the following table: 
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function 






interval 


monodromy 


interval 


monodromy 


T (i/ n )/ , 7 / t) 


2uii/n 


exp[-(2/n)?7i^ 


2a;i 


exp[— 2^im] 


2uj 


exp[-2^C] 


2nuij 


exp[— 2nr]jw] 


*? /B W) 


2ui/n 


exp[-(4/n)»7i£] 


u x 


exp[— 2t7xm] 


2ujj 


exp[-4?7^] 


nujj 


exp[— 2ra?7jM] 


4 1/n) M) 


2u 1 /n 


exp[-(6/n)?7i£] 


2wi/3 


exp[— 2^!^] 


2ujj 


exp[-6r/^] 


2nujj/3 


exp[— 2ra?7jM] 



Table 1: Monodromy Properties 
This means that, for example, 



x 



W»)(u + 2ui/n,{) = exp[-(2/n)r ?1 e]x( 1 ^(n,0, 



x 



(u,e + 2wi) = exp^Hx^KO, Vj = CM, J = l,2,3. 



The monodromies of the functions x^ 1//n - ) for n = 2,3,4 may be derived from the mon- 
odromy of the function x by using the following relations 

x{u^/2)x{u + u l ^/2) 



* (1/3) K0 



x^i/2) 

x(u, £/3) x(u + 2wi/3, f /3) x(u + 4a>i/3, f /3) 
x(2a; 1 /3,e/3)x(4u; 1 /3,e/3) : 



(5.18) 
(5.19) 



(i/4)/ n x(m, e/4) x(u + o;i/2, e/4) sfo + uu e/4) sfc + 3^/2, £/4) 
X WU " x(wi/2,e/4) i r(w li e/4)x(3w 1 /2,e/4) ' 1 " J 

The following factorisation identities of these solutions are useful (the dependence on £ 
on the l.h.s. is suppressed for brevity): 



x(u) x{—u) 
x^ /n \u) x {1/n \-u) 
4 ,n \u)x^ n \-u) 
x<?'*\u)x^\-u) 



. p (l/n) (u)+p (lA>) (e/n ; 



(5.21) 



, n = 2,3,4, 



(u) + p^(2^/n), n = 2,3,4, 



p^(«) + p wj,j (0 



(1/3), 



which give the potentials at various periods and are direct consequences of the first sum rule, 
as explained in ( |3.12| ). The second group is responsible for the renormalised constants in the 
extended twisted models (n = 1,2): 



x 



(!/»)/ 



(2u)x { j /2n \-u) + x^ n \-2u)x^ An \u) 



,(l/2n) , 



-^) (u) + p (l/n) (e/n) 
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-p^m/n), (5.22) 
x^ n \u)x^ 2n \-u) + x^ n \-u)x^ 2n \u) = -2p^ n \u) + p^ 2n \i/2n) 

+pW n \u 1 /n) 1 (5.23) 
4 /n \u)4 /2n \-u) + x d 1/n \-u)x d 1/2n \u) = -2p {1/n \u) + p^)(e/n) 

+p( 1 /")(^ 1 /n), (5.24) 

x(2u)x d m (-u)+x(-2u)x d m (u) = -pW(u) + pW(£/2). (5.25) 



They are obtained from the twisted second sum rules, ( |5.22| ) from ( |5.4j ), (|5.23[) from ( |5.5| 



and ( |5.24| ) from ( |5.6| ). Likewise we obtain the following relations from the twisted third sum 
rules, ( FTI2D and (|5TT3| ): 



x(3n)^ 1/3) (-m) + x(-3u)^ 1/3) (m) + x d (3u)a;J 1/3) (-2m) + x d (-3u)a;f /3) (2m) 

= -p^ 3 \u) + 3p(0-2p(2u Jl /3), (5.26) 

x(m)x (1/3) (-m) + x(-m)x {1/3) (m) + x(2w)xi 1/3) (-u) + x(-2w)4 1/3) (m) 

= -3p( M ) + p( 1 / 3 )(e/3) + 2p(2^ 1 /3). (5.27) 

5.3 Degenerate cases 

Twisted Calogero-Moser models can also be defined for the trigonometric and hyperbolic 
potentials. Since these potential functions have only one period, the half or third-period 
potentials are simply obtained as: 

2 ,, 2 f 4 cot 2 2m, 4 coth 2 2m 



cot u. coth u 



9 cot 2 3m, 9 coth 2 3w etc. 



The Lax pairs with spectral parameter for these potentials can be obtained from those Lax 
pairs given in section 4 by taking the degenerate limits of the functions {x, x^^^Xd, x d , ...}. 
With the help of the formulas (|5.18|) - (|5.20|) the degenerate limits of the twisted func- 
tions {a^ 1 / 2 ), x^/ 3 \ x d l ^ 2 \ x^ 3 \ ...} can be obtained from those of the untwisted functions 
{x, , Xd, x t , ...} given in section [O. 



5.3.1 Imaginary period infinite: Trigonometric potential 

In this case, as we have seen, the untwisted functions ( |3.27| ) 
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reduce to Q3.33Q 



x(u, £) = k (cot ku — cot k£) e a ^ u , Xd(u, £) = x(u, 2£), £) = x(w, 3£) k = Vl2a. 

From these we obtain, as expected, 

x (1/n) (u,0 = n/t(cotnra-cot<) e _< ", n = 2,3,4, 

xS 1/B) («,0 = x^\u,20, xl 1/3 \u,0 = x (1/3) (u,30. (5.28) 

In deriving these formulas from (|5.18|) - (|5.20|) use is made of an interesting trigonometric 
identity: 

!pi (cot(u + £) -cotu) , 

(cot w — cot v) I I — j; j-^ 5 r — = n (cot mi — cot ra;) , n = 1, 2, . . . , . (5.29) 

£=i l cot („) -cotv) 

The spectral parameter independent functions are 

x(u) = Xd{u) — xt{u) — a (cot au =F 1), = '(u) = 2a (cot 2au =p 1), 

X (V3) ( M ) = X (V3) (m) = x (l/3) (m) = 3a (cQt 3qu T !) ; 

x (1/4 )(m) = 4 1/4) (w) = 4a(cot4cra^l), a : const. (5.30) 
5.3.2 Real period infinite: Hyperbolic potential 

In this degenerate limit, twisting along the direction {2u>i, 2u> 3 } — > {2ui/n, 2u> 3 } gives trivial 
results: 

x (u) = x {1/2 \u) =x (1/3) (m), x d (u) =xi 1/2) (w),..., (5.31) 

since c^i = oo ( |3.25| ). That is, the twisted models are identical with the untwisted models. 
Twisting along other directions give non-trivial results. For example, 

{2u;i,2u; 3 } -> {2u u 2u 3 /n}, 

the degenerate limit of the functions can be obtained in a similar way as above: 

x(u, £) = k (coth ku — coth k£) e^, xa(u, £) = x(u, 2£), £) = x(w, 3£) k = vl2a. 

From these we obtain, as above, 



x (1/n) (u,0 = (coth riKu- coth «£) e a5 ", n = 2,3,4, 

4 Vn) K0 = ^)( M ,20, xf /3) K0 = x (1/3) K30- (5.32) 
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The spectral parameter independent functions are 

x(u) = x d (u) = x t (u) = a (coth.au ± 1), x^ 2 \u) = \u) = 2a (coth2a-u ± 1), 
x<V 3 \u) = x d 1/3) (u) = 4 m (u) = 3a (coth3au ± 1), 

x [1/A) {u) = x d m (u) = 4a (coth4au ±1), a : const. (5.33) 

6 Summary and comments 

The construction of universal Lax pairs for all of the Calogero-Moser models based on root 
systems is completed. This paper presents universal Lax pairs for the twisted models based 
on non-simply laced root systems, following those for the models based on simply laced root 
systems given in the first paper [|IJ and those for the untwisted models based on non-simply 
laced root systems constructed in the second paper ||. All of the Lax pairs for the twisted 
models presented here are new, except for the one for the _F 4 model based on the short roots. 
These Lax pairs are for the elliptic potentials and contain a spectral parameter. The explicit 
forms and the properties of the functions appearing in the Lax pairs, the untwisted and 
twisted functions, are determined from the functional equations (which we call sum rules) 
necessary and sufficient for the consistency of the Lax pairs. 

The spectral parameter dependent Lax pairs for the hyperbolic, trigonometric and ratio- 
nal potentials for the untwisted as well as the twisted models are derived as the degenerate 
limits of the elliptic potential case. These, we believe, have not been reported yet. All of 
the Lax pairs discussed in this paper are of the root type and contain at least as many 
independent coupling constants as the number of independent Weyl orbits in the root sys- 
tems. Our main result is very simple. The Lax pairs for the twisted models can be obtained 
from those for the untwisted models (given in paper II) by replacing the untwisted functions 
{x,y,z, x d ,y d ,z d , x t ,y t ,z t } by the twisted functions {x (1/n) , y (1/n) , z (1/n) , x d /n \y d /n \ z d /n \ 
Xf 1 , yt i z^ n ^}. This also applies to the other type of universal Lax pairs. The minimal 
type Lax pairs can be obtained simply from those for the untwisted models by replacing the 
untwisted functions {x, y, z} by the twisted functions {x^^ 2 \ y^^ 2 \ z^ 1 ^} (for the Gi model, 
{a^ 1 / 3 ), y^ 3 \ z^l^Y) for the short roots. 

As for the twisted models based on B n , C n and BC n root systems, a new type of potential 
terms with independent coupling constants can be added without destroying integrability. 
They are called extended twisted models and have three, three and five independent coupling 
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constants for B n , C n and BC n models, respectively. The Lax pairs for the twisted G2 model 
have new features. This model has an elliptic potential function with a pair of primitive 
periods {2u\, 2co> 3 } for the long roots and another elliptic potential function with third-period 
{2a>i/3, 2uj^} for the short roots. New sum rules appear in the proof of the consistency of 
the Lax pairs for this model. 
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